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Abstract. The Dirac operator enters into zero curvature representation for the 
cubic nonlinear Schrodinger equation. We introduce and study a conformal map 
from the upper half-plane of the spectral parameter of the Dirac operator into 
itself. The action variables turn out to be limiting boundary values of the imagi- 
nary part of this map. We describe the image of the momentum map (convexity 
theorem) in the simplest case of a potential from the Schwartz class. We apply 
this description to the invariant manifolds for the nonlinear Schrodinger equation. 



1. Introduction. The nonlinear Schrodinger equation (NLS) with repulsive 
potential* 

W = -^" + 2MV, 

has been solved for square integrable initial data on the circle and for a sumable 
potential on the entire line, see e.g. [MCV, FT]. It is an important problem 
which comes from statistical mechanics [MCV] to prove complete integrability 
for extensive class of initial data which accommodate some translation- invariant 
Gibbs' distribution on the entire line. 

The integration of the NLS equation on the entire line is based on the scat- 
tering theory for the Dirac operator — the simplest 2x2 canonical system. 
In this paper we introduce a special conformal map (CM) from the upper-half 
plane of the spectral parameter into itself. The action variables, known before 
for the rapidly decreasing case in terms of so-called scattering coefficients, are 
defined now as a boundary values of the imaginary part of this map. The limiting 
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boundary values are not a problem: they are known always to exist for a positive 
harmonic function [KO]. What we propose, in fact, is a continuum analog of the 
"Marchenko-Ostrovskii" conformal map [MO]. 

It is known that the Dirac operator is self-adjoint and that the Weyl- 
Titchmarsh functions are unique under sole condition of continuity of the po- 
tential. We establish a relation between the boundary values of the CM with 
the Weyl-Titchmarsh functions, with the hope that, suitably interpreted, this 
relation defines a CM for sufficiently broad class of potentials. 

Construction of actions for the NLS equation with initial data from a par- 
ticular function class reduces to the description of all possible boundary values 
of the CM, i.e. to the image of the momentum map. In this paper, we prove 
the simplest theorem of that type: for the Schwartz' class S(R 1 ) of rapidly 
decreasing infinitely smooth complex functions, the image of the momentum 
map is a convex cone of real non-negative functions likewise from S^R 1 ). This 
is an infinite-dimensional generalization of the convexity theorem of Atiyah- 
Guillemin- Sternberg [AT, GS]. As an application of our result we solve the fol- 
lowing problem. 

Problem. What values of conventional local integrals of motion H±, H2, • • • , 
if TV, N < 00, where 



correspond some invariant manifold in S(R 1 ) ? What are all the invariant man- 
ifolds in S(R 1 ) with given values of Hi, H2, ■ ■ ■ , Hjy ■ 

The paper is organized as follows. In sections 2 and 3 we present standard 
facts concerning Hamiltonian formalism and scattering theory needed in our 
considerations. The conformal map itself is introduced in section 4. It turns out 
conventional integrals of NLS are the moments of the positive measure defined 
by the imaginary part of the CM. Section 5 establishes a relation between the 
boundary values of the CM, Weyl-Titchmarsh functions and the spectral matrix. 
Section 6 contains necessary information about the periodic problem for the NLS 
equation. The relation between the CM and "Marchenko-Ostrovskii picture" is 
presented in section 7. We also consider the whole structure in the linear limit. 
An interpretation of the results from the point of view of scattering theory is 
given in section 8. In the two sections 9 and 10, we show that the actions are 
canonically paired to the angles, which are given in terms of the phase of the 
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reflection coefficient. This demonstrates an action of the torus with continuum 
number of generators, a standard requirement in the convexity theorems. We 
exploit this structure to give an effective solution of the moment problem. The 
convexity theorem is presented in section 11. A conjecture about the image of 
the momentum map for potentials from the Sobolev' class H s , s > is also 
stated there. The last section 12 contains the solution of the problem. 

The author is grateful to P. Deift, L. Faybusovich, H. McKean, V. Peller and 
N. Young for discussions. 

2. Hamiltonian formalism and zero curvature representation. * We 

denote by S^R 1 ) the Schwartz' space of complex rapidly decreasing infinitely 
differentiable functions on the line such that 

sup|(l + x 2 ) n / (m) (x)| < oo forall n,m = 0,1, ■■■ . 

X 

On the phase space M = S(R 1 ) x S(R 1 ), which is a product of two copies of 
Schwartz' space, we introduce a Poisson bracket 



{A, B}=i [ 



OA dB OA OB , 
— = ——dx. 



R i dip dip dip dip 

The bracket is well-defined for any functionals A, B with "nice" gradients. The 
original NLS equation 

iip* = -iP" + 2\iP\ 2 iP, 

can be written in the form 

ip" = {8H 3 , iP} with 8H 3 = / l^'l 2 + \ip\ 4 dx = energy. 

Jr 1 

There are two other classical integrals 2Hi = J \ip\ dx =number of particles and 
4H2 = — i J ipip'dx =momentum. Hi, H2 and H% are the first three Hamiltoni- 
ans of an infinite hierarchy of commuting flows. NLS is a consistency condition 
for the zero-curvature representation 

where U = Uq + \U\ with 

V. = I I fj , U 1 = -ia 3 



*The commutator representation for NLS equation was found by [ZS] ; see also [AKNS] for 
a more systematic approach. 
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and V3 = Vq + AVi + \ 2 V 2 with 



Vo = i( l f t. 2 )> V 1 = -2U 2 , V 2 = -AU,. 



Here and below a denotes the Pauli matrixes 



o n / -A MO 



;Tl > i o J • a -' ^ i o ) ■ a:i v<» 1 

3. The reduced monodromy matrix. The results of this section are basically 
well known. They are included here to make the paper self-contained.* The 
standard monodromy matrix M(x, y, A) is the 2x2 solution of 

M'(x,y,X) = U(x,\)M(x,y,\), M(x, y, X)\ x=y = /, 

where / is the identity matrix. The reduced monodromy matrix T(x, y, A) is 
defined as 

T(x, y, A) = E-\Xx)M(x, y, X)E-\-Xy), 

where E{Xx) = exp(— iXxa^) is a solution of the free equation (ip = 0). The 
matrix T solves 

T'(x,y,\) = U (x) E(2Xx)T(x,y,X), T(x, y, X)\ x=y = I. 

It can be expressed as 

T(a;,y,A)=exp / U ^)E(2X^)d^. 
•Jy 

Note that both M and T belong to SX(2,C). The symmetry of the matrix 
U : aiUo(x)o~i = Uq(x) is inherited by T : a\T(x, y, A)oi = T(x, y, A). 
For ip E S(R 1 ) and A real, the existence of the limit 

limT(x, y, A) = T(A) = (jfy ^ , for y -> -oo and x -> +oo 

follows easily from the inequality 

|exp f U (Z)E(2\Z)dt-I\ <exp ! \U Q {£)\<% - 1, 



*A study of the scattering theory for canonical systems goes back to [K,KMA,MA]. For the 
latest results see [BDZ]. 
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where | • | is any multiplicative matrix norm. Due to the symmetry of T(x, y, A), 
T(A) has the form 



The important fact needed bellow is that 6(A) G S(R 1 ). This implies the follow- 
ing property of a(A) 

i. For A real |a(A)| > 1 and |a(A)| 2 - 1 G S{R X ). 

To establish properties of T(A) for complex A, we introduce the Jost solutions 
M±(x,\): 

M' ± (x, A) = U(x, X)M±(x, A), M±{x) = E(Xx) + o(l) as x -> ±oo. 

They can be expressed through the solution of the free equation by means of 
transformation operators: 



Integral equations and estimates for the kernels F± is 2 x 2 are derived following 
[M, section 3.1]. 

The kernels r± are unique and infinitely smooth in both variables. They have 
a symmetry T± = aiT±a\ and admit unique decompositions T± = C± + A±, 
where C±, resp. A± commute, resp. anti-commute with o\. The components 
C_, A- of the kernel T_ satisfy 






X 





— oo 



x 




t+x 
2 



whence 



X 





— oo 
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The components C + , A + of the kernel T + also satisfy 

+00 

C+(x,t) = - J U (Z)A + (Z,t + Z-x)dt, 

X 

t + x 

A+(x,t) = -^U ( t -^)- J U (t)C+(t,t + x-t)dt, 

x 

whence 

A+(x,x) = ~U (x), C + {x,x)= l - J U*(£)dt. 

x 

The monodromy matrix M(x, y, A) can be written in the form 

M(x,y) = M + {x)M- 1 {y) = M_(x)MZ 1 (y). 

Therefore 

M+ 1 (x)M(x,y)M-(y) = M- 1 (y)M-(y) = M+ x (x)M-(x). 

The expression at the right does not depend on x or y and it is equal to T(X) 
as one can see by passing to the limit as x — > +00, y — > —00. Introducing the 
notation M± = (Mj_, M|) for the columns, we get a(A) = det (Mi, M\). Now 
the following properties of a(A) follow easily from the integral representation of 
the Jost solutions 

ii. a(A) is analytic in the upper half-plane and infinitely smooth up to the 
boundary; 

in. a(A) = 1 + o(l) as |A| — > 00 and a(A) is root-free. 

4. The conformal map Poo(A). Let Poo(A) be such that a(A) = exp(— zpoo(A)) 
for A in the upper half-plane. From the properties of a(A) presented above follow 

i. Poo (A) is analytic in the upper half-plane and infinitely smooth up to the 

boundary. 
ii. %oo(A) > for 3A > 0. 

This means simply that Poo(A) defines a conformal map from the upper half-plane 
of the spectral parameter into itself. Moreover, 

Hi. Poo(A) = o(l) for |A| — > 00; for real A, the density of the measure 
^oo(A) = Qpoo(A) dX belongs to S^-R 1 ) and dfx 00 (X) has moments of 
any order. 
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This is due to the fact that |a(A)| 2 — 1 G S(R 1 ). The function £>oo(A) can be 
written in the form 

1 f dfiooit) 

Poo (A) 



TV I t - A 



see [A]. Expanding the denominator in inverse powers of A, we obtain: 



°° 1 1 f +oc 
Poc(A) = -^- TT - / t k d^(t) = 



Hi H2 H3 
T~ A^~ A 3 " 



where Hi,H 2 and #3 are the integrals introduced above. The expansion has 
an asymptotic character for A : 5 < argA < tv — 5. An identification of the 
coefficients of the asymptotic expansion employs a recurrence relation for H's. 
As in [MCV], a straightforward computation produces: 



2iA ^ , ; ; + D- 



dip(x) dip(x) 



, dT -dT 

dip dip 



2ip 



2zA^ 



dip(x) dip(x) 
Here we have used only the formulas for derivatives: 



, dT -dT 

dip dip 



dT 

dip(x) 



= T(+oo 



0) 



E(2Xx)T(x,-oo), 



^=T(4co,»)(JJ J)B(2A,)r(x,-oc). 
Now we can put the formulas into desired form of a recurrence relation: 



JVH n+1 =KVH n , n = l,2, 



with" 1 



K = 



-iipD~ x ip -i /2D + iipD~ x ip 



-i/2D + iipD~ x ip 



-iipD ip 



and J = ioi- 



This allows us to reconstruct the if's starting from 2H\ = J \ip\ . The same 
recurrence relation implies an involutive character of i/'s; see [MCI]. 



*d- i = i/2 (/foo-zr)- 
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5. The relation between Poo(A), Weyl-Titchmarsh functions and the 

hi 



spectral matrix. By definition, the Weyl solutions h± = ^ J of h! = Uh 

belong, resp., to L 2 (R±) for A with strictly positive imaginary part; they are 
normalized by h+(x, A) = 1, resp., h]_(x, A) = 1 for x = 0. They can be written 
in the form 

h+(x, A) = M 1 (x, A)m+(A) + M 2 (x, A), 
h-(x,X) = M 1 (x, A) + M 2 (x, A)m_(A), 

where M(x, A) = M(x, 0, A). For example, for the free equation, the Weyl- 
Titchmarsh functions m± are identically zero and 



h+(x,X)= (JL) £L 2 (R\), 



g i\x 



h-(x,\)=\^ Q j eL 2 {R\_). 

It is well-know fact, see [LS], that under the sole assumption of the continuity 
of the potential tp, the functions m± are unique. 
Our first goal is to prove that 



e 





jl - 


m_m+ 


2 


(1-1 


m_|_ 2 )(1 — 


m_ 


2 ) 



29p - (A) = , ™, + ^ for A real. 



This formula provides an alternative way to define the CM. 

Note that M 2 resp., M\ are from L 2 (i?^_) resp., L 2 (i?^_) for A with strictly 
positive imaginary part. Therefore 

M\ = k-(X)h-, Ml = k+(X)h+, 

for some numbers k±(X). Take ip supported in [—L,+L]. Then 

M(0,L)Ml(L) = AfJ(O) = k+h-(0). 



In more detail, 



z :£H(A)=(7K 



or 



m22e lXL = k+, mi2e lXL = m+k + . 
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Finally, for A real 

1 = \rri22\ 2 — |^i2| 2 = |^+| 2 (1 — l m +| 2 ) an d \k+\ 2 = 
and likewise 



1 — |m+| 



\k-\ 2 = 



1 



1 — |m_ | 2 
Therefore, for such A, 

e 29 Poo (A) = | a ^^|2 = | det ( M l ^ M 2^|2 = j^^^ _ TO _ m+ ) 



11 — 171-171- 1 J 



(1 — |m_|_| 2 )(l — |m_ | 2 ) ' 

as required. 

The first order system i 7 " = t/F introduced at the beginning can be expressed 
as an eigenvalue problem for the Dirac operator £): 



lj + l# 



/ = A/. 



The Spectral Theorem. For t/ie operator D with a continuous potential there 
exists a unique 2x2 Hermitian spectral matrix ||py(A)||; \\pij(X') — Pij(\")\\ is 
positive definite for any real A' > A"; the entries Pij(X) are functions continuous 
from the right and 

\dpij(X)\ 



1 + A 2 



< oo. 



For any vector-function f E L (R ), the sequence 



f T (x)M (x,X)dx, i = 1,2 

-n 

converges to some Fi(X) in L 2 (i? 1 , pij): 

f " FKFi - F')d P y 

J ij 

Moreover, 



oo, as n — > oo. 



\f\ 2 dx= j £F*(A)Fi(A)dp y (A). 



10 



K.L. VANINSKY 



The function pa, i = 1,2 is the boundary value of the expression 



1 1 




m_ 


2 


m_|_ 


2 


2tt 


1 - 


- m. 


-TTl + 


2 



and pi2 — P21 is the boundary value of the expression 



1 m+(l - 


m_ 


2 )-m_(l- 




2 ) 


2tt 


1 - 


- m_m+ 


2 



1 (1-1 


m + | 2 )(l-| 


m_ 


2 ) 


4tt 2 


1 — m_m + 


2 



The proof is by standard arguments, see for example [CL] . By direct compu- 
tation, 

detHpiJ 

which implies: det||py(A)|| = e ~ 2Qp °°^ /4tt 2 for A real.* 

One can think about spectral matrix as a smooth functional on the phase 
space. Indeed, gradients of the functions m± can be computed, e.g. 

We split the proof of the first formula into 5 steps. 
Step 1. By direct computation 

W °Wo), 



dip(y) 
where < y < x. 



Step 2. Let /, g be two vector solutions of the equation: Df = A/, then 
det[f(x),g(x)] does not depend on x. The proof is trivial. 

Step 3. The Weyl solution h + normalized by: h+(x,y, \)\ x=y = 1, can be 
expressed 

h+(x, y, A) = M\x, y, A)m+(A, y) + M 2 (x, y, A), 



*H.P. McKean in [MC2] proved a similar formula. 
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where m + (X, y) is the corresponding Weyl-Titchmarsh function. In particular, 
h+(x, A) = h+(x, 0, A). The fact needed below is that 

fr+(y, 0, A) 

To prove this, note that the identity M(x, 0) = M(x, y)M(y, 0) implies 

M^x, 0) = mn(y, OjM^x, y) + m 2 i(y, 0)M 2 (a;, y), 
M 2 (x, 0) = mi 2 (y, 0)M 1 (x, y) + m 22 (y, 0)M 2 (x, y). 

Therefore, for < y < x: 

h + (x, 0, A) = M 1 (x, 0)m+(A, 0) + M 2 (x, 0) 

= [m n (y, OjM^a;, y) + m 21 (y, 0)M 2 (x, y)] m+(A, 0) 

+ [m 12 (y, OjM 1 ^, y) + m 22 (y, 0)M 2 (;r, y)] 
= [m 2 i(y, 0)m + (A, 0) + m 22 (y, 0)] 



x 



(X,?/j m 21 (y,0)m + (A,0) + m 22 (y,0) +M ^ 



The result follows. 

Step 4- The purpose of this step is to prove the formula 



<9V>(y) 
where 



dm + (A,0) 

= -m + (A, 0)m+(\,y) 



mi2(y,0) mn(y,0) 
A B 



A= m 12 {y,0) - m + (A,y)m 22 (y,0), 
5 = -mn(y, 0) + m+(A, y)m 2i (y, 0). 

Consider the eigenvalue problem: Df = Xf with the boundary conditions 
/ 2 (y, A) = 1, / 1 (6, A) = f 2 (b, A), where < y < 6 < +oo. The solution can be 
written in the form 

f(x, A) = M 1 (x, y, A)m+(A, y, b) + M 2 (x, y, A). 

Due to the selfadjointness of the spectral problem, the limit 

m+(X,y) = lim m+(X,y,b) 
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exists and does not depend on the boundary condition at the right end of the 
interval. From the definition, 



Therefore, 



m+{X,y,b) = (o, y, A). 

mn - m 2 i 



Vm + (A,0,6)= m22 - mi2Vm22 ' Vmi2 (6,0,A) 



mn - m 2 i m 22 - m 12 
m 22 - mi2 Vmn - Vm 2 i 
mn - m 2 i mn - m 2 i 



(6,0, A). 



By Step 1, 



dm+(A,0,6) 



m+(A,0,6) — — 

(m 22 - m 12 ){b, 0) 

f\ n t\ (^12 m 22 ){b : y)mn(y, 0) 

m+(A, 0, 6) — — 

(mn -m 2 i)(o, 0) 

m+(A, 0, 6)(m 2 2 - mi 2 )(6, y) 

mi 2 (y, 0) mn(y,0) 



x 



(m 2 2 - mi 2 ) (6,0) (mn - m 2 i)(6, 0) 
Using the identity M(6, 0) = M{b, y)M(y, 0), and simple algebra, one finds 



-m+(A,0,6)m+(A,y,6) 

"7-12 (y,0) 



+ 



mn(y, 0) 



.mi 2 (y,0) -m + (X,y,b)m 22 (y, 0) -m n (y, 0) + m+(A, y, 6)m 2 i(y, 0) 

Now pass to the limit 6 — > oo. 
5iep 5. By Step 3, 

m + (A,y)[m 2 i(y,0)m + (A, 0) + m 22 (y, 0)] = mu(y, 0)m+(A, 0) + mi 2 (y,0). 
After simple algebra, 



m + (A,0) _ 1 
A ~ B 



Therefore, 



<9m+(A,0) 



= -m+(X,y) 



m + (A, 0)m n (y, 0) + m 12 (y, 0) 
-mu(y, 0) + m 2 i(y, 0)m+(A, y) 
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By Step 3, 

<9m+(A,0) h\(y,X) 



mil(2/ ' 0) /^o)- m2l(2/ ' 0) 



mu(y, 0)^(y, 0) - m 21 (y, 0)^(y, 0) ' 

By Step 2, the denominator can be computed at y = 0, where it is equal to 1. 
We are done. 

The proof of the second formula is the same. 

6. The necessary information about periodic problem*. In the next two 
sections we relate the CM defined by Poo(A) with the analogous object known in 
the periodic case. Consider the periodic problem for the NLS, i.e. x is on the 
circle of perimeter 2L and assume that the function ip is infinitely smooth. The 
periodic/ anti-periodic, spectrum of the auxiliary spectral problem is determined 
by 

Ax, (A) = traceM(L, —L, A)/2. 

Namely, Ao is a periodic, resp. anti-periodic, eigenvalue if it is a root of the 
equation: A L (X ) = +1, resp. —1. The function 7t L (X) is defined by: A L (X) = 
coshz7Ti / (A); it has the properties 

i. tvl(X) is single- valued on the A-plane cut along open gaps (A 6 E 1 : 
|Al(A)|>1); 

ii. tv l (X) = tt l (A); 

Hi. ^s7Tl(X) > for A in the upper half-plane. 

The function k~l(X) determines a conformal map from the upper half-plane into 
a slit-domain. In fact, this CM was known to physicists already in 1929, see [B]. 
This map was used by [MO] in their description of all possible spectral sequences 
for the periodic Schrodinger operator. 

The function ttl(X) can be written in the form 



tt l (A) = 2XL + p L (X) 



in which 

i. Pl (A) is analytic and continuous up to the boundary; 
ii. pl (A) =pl(A); 

Hi. ^spl(X) > for A in the upper half-plane; the inequality is strict, except 
for a simple zeroes, one in each band. 



*For more information see [MCV]. 
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Introducing a measure <Ihl{\) = $spl(\)cI\, we can write £>l(A) in terms of the 
Cauchy integral 

t-X 

The function Pl(A) has an asymptotic expansion at infinity analogous to the 
asymptotic expansion of Poo(X): 

°° 1 1 p + OO TT TT TT 

MA) = . l( A)-2AL = -i: / ^(0 = -^-^-^-.... 

fc=0 >/_00 

7. The infinite-volume limit. Deformation to the linear problem. The 

scaling ip — > v / xV'? H > 0, transforms the original NLS equation into 

fi/,* = -^" + 2x|^| 2 <i/>- 

We denote by p^(A), ^(A), etc., all the quantities introduced before to em- 
phasize their dependence on the parameter x. The function p^(A) turns out to 
be a continuum analog of p£(A). Namely, for A with strictly positive imaginary 
part 

Pl(A)^p£(A), asL^oo. 

To prove this fix a compactly supported potential. Then T(A) is an entire matrix- 
function, as follows from the matrix exponential. Consider now the identity 

T(L, -L, A) = E(XL)T(X)E(XL) = ( „*f£j*XL 



6(A) a*(A)e^ 

Taking the trace of both sides 

2A£(A) = e ^ 2XL+ ^^ + e -^+^W) = e - 2lXL a(X) + a*(X)e 2lXL . 

Therefore, the difference e'^W - e -iplW = - e 4iAI, [e^W-|-a*(A)] is negligible 
as L — > 00. That gives us a second proof that Poo(A) has a positive imaginary 
part. 

The deformation of the measure d/i* in the infinite-volume or linear limit can 
best be expressed by a " commutative" diagram: 



^- 1 dMS,(A),^|0 



n 

where 

$(n) = \ ! + e- ix ™/ L iP(x)dx 7 V>(A) = / e~ ixX ^{x)dx. 
The convergence of measures has to be understood in the weak* sense. 
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8. The scattering theory. This section is a brief retreat into physics. It 
explains on physical grounds why Poo(X) has a positive imaginary part. Consider 
two solutions f\ and f r with prescribed asymptotics as x — > ±00: 



where 



x — 00 + 
/, /_ +J B(A)/_ A(A)/U 



< iXx \ ( 



/- = ( 1 ' ; - 



are the solutions of the free equation. The solutions //, resp. f r , describe a 
particle which approaches a potential barrier from the left, resp. right. The 
coefficients B, resp. C, are called left, resp. right reflection coefficients, A = D 
is a transmission coefficient. 

As in sect. 5, we can compute a scattering matrix S'(A): 



5(A) 



A B\ ( l/a(A) -6(A)/a(A)\ 
C D ) ~ \b{\)/a{\) l/a(A) J 



Compute probability for the particle to go through or been reflected by the 
barrier: 

Prob {particle goes through the barrier} = / \A(X)\ 2 \(f(X)\ 2 d\ 

Prob {particle is reflected by the barrier} = / |S(A)| 2 |(/(A)| 2 rfA, 

Jr 1 

where /(A), f |/(A)| 2 rfA = 1, is an arbitrary function, see [LL]. Since Prob{-} < 
1 with necessity we obtain |A(A)| = |e iPoo( ^ A ^| < 1, for all A and therefore 
3poo(A) > 0. That is a third proof that Poo(X) has a positive imaginary part. 

9. Action of H's flows on 6(A) and the angles. Now we return to the main 
topic. In the next two sections we demonstrate that /(A) = Qpoo(A) produces on 
invariant manifold a torus action with continuum number of generators. Note 
that the action-angle variables for NLS equation in the scattering case were 
constructed in [ZM]. We prefer here a method of computation which uses the 
CM and is therefore different from the conventional one. 
By direct computation, we obtain 



b'{X) = {H 1 ,b(X)} = -i/2b{X). 
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Another useful fact is 



KVb(X) = AJV6(A). 



Then, the inductive rule 



{H n+1 ,b} 




implies 



6- (A) 



{H n+1 ,b} 



i/2X n b(X). 



This means that 6(X) = —phase 6(A) changes with a uniform speed; it is an 
angle on the invariant torus. 

10. Effective solution of the moment problem. To show that /(A) is 
canonically conjugate to 0(A), we establish a formula for an effective solution 
of the moment problem. This was introduced in [V] under some restrictive 
assumptions. Here we relax the assumptions and present an easy proof, which 
can be generalized to the continuum case.* 

The discrete case. Let H n be a sequence of moments of some positive measure 
Hk on Z 1 : 



converges absolutely for all k, then the measure is unique and nu = Sk for all k. 
It is easy to prove the formula for compactly supported [i. Take the function 

sin7r(x — k) 



*For another approach to effective solution of the moment problem on the finite interval 
see [F, chapter 9, volume 2]. 




k 



If the series 




tv(x — k) 
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which vanishes everywhere on Z 1 



except x = k and consider its Taylor expansion 



sin7r(x — k) 



oc 



x sin7r(x — k) 



l(n) 



tv(x — k) 



E 




n=0 



x=0 



Integration of both sides with respect to \i produces the formula. To prove it 
under the stated assumptions one has to consider the even and odd parts of the 
measure separately: e n = (fx n + / u_ n )/2, o n = (fi n — (i- n )/2. Take the even 
part. All moments of odd order vanish and the even order moments match with 
the corresponding moments of the original measure. The formula is true for the 
measure with a cut-off, which now can be removed by the bounded convergence 
theorem. The odd part can be considered in a similar way. The final result is 
produced by taking a sum of both parts. The proof is finished. 
Denote by (5 a the Gaussian density 



and by f a = f <E> 17 the convolution of a function or measure / with the kernel 



The continuum case. Let H n be a sequence of moments of some measure on 
R 1 . If the series 



converges absolutely for all x and some a > 0, then the measure is unique and 
s a (x) = n a (x) for all x. 

The proof is identical to the discrete case, only now one has to integrate the 
Taylor expansion of — y). 

Introducing /^(A) and using the preceding formula 



(5 a . 




n=0 



{2I°(\'), 6>(A")} = ^ - A )tlo {2#»+i, W)} 



°° (\"\ n , ^ 

= J2 A l-^ a ( X> -A)]f;l = ^(A' -A"). 

n=0 U - 



Passing to the limit as a | 0, we obtain 



{i(\'),e(\")} = 5(\'-\"). 



The involutive character of the Ps can be proved in a similar way. 
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11. The Convexity Theorem. For a potential from the Schwartz' class 
S(R 1 ) of complex rapidly decreasing infinitely smooth functions, the image of 
the momentum map is the convex cone of all nonnegative functions from the real 
part of S(R 1 ) . 

The part that the map is into was already explained in section 4. The onto part 
is more difficult. Assume that 6(A) is from ^(-R 1 ) and pure real for definiteness. 
Then the theorem is proved using standard arguments, by reconstructing the 
potential from the scattering matrix, see [M]. This involves a solution of the 
Gel'fand-Levitan integral equations for the kernels T± of the transformation 
operators: 



T + (x, t) + tt(x + t) + / T+(x, s) Q(s + t)ds = 0, 



where 



+oo 



T-(x,t) + n(x + t) + / T-(x,s)n(s + t)ds = 0, 



and 



where 

_ +oo _ 

fi <*>=U) 5( o X) )- *(*> = s/™^ 'W = "S- 

— oo 

The infinite smoothness and decay of the kernels T± are provided by the fact 
that Q are from S(R 1 ). 

The theorem implies that the conditions in section 4 which describe the func- 
tion Poo (A) for potentials from the class S^-R 1 ) are not only necessary, but also 
sufficient for Poo(A) to be produced by a potential from S(R 1 ). 

We conjecture, by analogy with the Fourier transform (see section 7) , that for 
a potential ip from the Sobolev class H s , s > 0, i.e., (1 - A)~ s / 2 tfj e L 2 (i? 1 ), the 
image of the momentum map fills up the convex cone of nonnegative functions 
of the weighted space L 1 (i? 1 , (1 + A 2 ) s dX). 

12. Solution of the Problem. The only condition on a finite number of 
Hi , • • • , Hn to correspond some invariant manifold is 
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In fact, it is the condition for a truncated moment problem to be definite, see 
[A]. 

A description of the invariant manifolds in S(R 1 ) with given values of Hi, ■ ■ ■ , 
Hn reduces to a construction of the function Poo(X) which satisfies the conditions 
i.-iii. of section 4 and has the prescribed asymptotic expansion 



Poo (A) = - 



Hi 
A 



#2 

A 2 



Hn _ / 1 s 
'W~° [ X^ ) 



for 5 < arg A < tt — 5, 

We will give a description of all such functions with a given value of just one 
integral Hi : 

where the arbitrary function y?(A) is such that 

%. if(X) = o(l) for |A| — > oo within the angle 5 < arg A < tv — 5; 
ii. Scp(X) > for 3A > 0; 

Hi. Q(p(X)/\X + (p(X)\ 2 is of the class S(R 1 ) on the real line. 

Indeed, according to [A, section 3.4], any function Poo(A), such that Qpoo(A) > 
for SA > and with the asymptotic behavior 

Poo(A) = — ^ ~o(j) 

within the angle 5 < arg A < rr — 5 can be represented in such a form with an 
arbitrary <f(X), satisfying i.-ii. and vice versa. The proof is easy. Consider 



<p(\) = -\- 



H 



Poo (A) ' 

This implies the condition i. Now 

S<p(\) = -3A + 3A J 

By Schwartz inequality 



With Poo (X) = 1 J J dfMooit) = TV Hi . 



dfx c 



ttH, 



t-X\ 2 



J t-X 



2 • 



t-X 



< 



djJL c 



dfx c 



It- A| 



txHi 



It- A| 



and this implies the condition ii. The proof of converse statement is trivial. In 
order to get Poo(A) which corresponds to some invariant manifold in S(R 1 ) we 
need to impose some restrictions on (p(X). For real A, 



^Poo(A) 



ffiSy(A) 
\X + <p(X)\< 
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This produces the additional condition Hi. Obviously, there are plenty of suitable 
functions <p. 

The case of general N < oo can be considered in a similar way with an aid of 
the so-called Schur algorithm, see [A, section 3.3]. 
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